GW quasiparticle band structures of stibnite, antimonselite, 
bismuthinite, and guanajuatite 
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We present first-principles calculations of the quasiparticle band structures of four isostructural 
semiconducting metal chalcogenides A2B3 (with A = Sb, Bi and B = S, Se) of the stibnite family 
within the Go Wo approach. We perform extensive convergence tests and identify a sensitivity 
of the quasiparticle corrections to the structural parameters and to the semicore d electrons. Our 
calculations indicate that all four chalcogenides exhibit direct band gaps, if we exclude some indirect 
transitions marginally below the direct gap. Our converged quasiparticle band gaps are 1.6 eV 
(Sb2S3), 1.3 eV (Sb2Se3), 1.7 eV (Bi2Ss) and 1.2 eV (Bi2Se3). By comparing our calculated gaps 
with the ideal Shockley-Queisser value we find that all four chalcogenides are promising as light 
sensitizers for nanostructured photovoltaics. 

PACS numbers: 71.20.-b 74.70.Xa 78.20.-e 78.56.-a 



I. INTRODUCTION 

The development of sustainable energy solutions based 
on scalable processes and non-toxic materials constitutes 
a key priority in the current scientific research agenda, 
and in this area nanostructured energy-harvesting solar 
and thermoelectric devices are playing a lead role. Re- 
cently there has been a surge of interest in devices using 
semiconducting metal chalcogenides of the stibnite fam- 
ily. For example recent studies have demonstrated the 
potential of these semiconductors both in photovoltaics 
applications,— ~— and in thermoelectric generators.— 

In the area of nanostructured photovoltaics semiconduct- 
ing metal chalcogenides have successfully been used to 
replace the inorganic dye in dye-sensitized solar cells^ 
leading to the development of solid-state semiconductor- 
sensitized solar cells.—— In these devices thin layers or 
nanoparticles of the semiconducting chalcogenides act as 
light absorbers, and upon photoexcitation they transfer 
an electron to the acceptor (typically Ti02) and a hole 
to the hole-transporter (for example a conducting poly- 
mer). The record efficiency within this class of devices 
is 5.1% and was obtained using stibnite (Sb 2 S3) as semi- 
conductor sensitizer X 

A recent atomistic computational study of photovoltaic 
interfaces for semicondictor-sensitized solar cells pointed 
out that, in addition to stibnite, the other members of the 
stibnite mineral family, namely antimonselite (Sb2Se3), 
bismuthinite (B12S3), and guanajuatite (Bi 2 Se3), exhibit 
optical properties similar to stibnite and should be con- 
sidered as potential candidates for novel semiconduc- 
tor sensitizers £ Using density-functional calculations and 
empirical scissor corrections of the band gaps, in Ref|H it 
was found that stibnite and antimonselite should form 
type-II heterojunctions with Ti02, while bismuthinite 
should form a type-I heterojunction and hence would not 
be able to transfer electrons to Ti02- These theoretical 
predictions have recently been confirmed by the exper- 
imental investigations of RefsJJH, thereby providing a 
motivation for further studies and for the more sophisti- 



cated analysis presented in this work. 
The four minerals of the stibnite family crystallize in 
an orthorhombic structure consisting of parallel one- 
dimensional (A4Be)„ ribbons, with A = Sb, Bi and B 
= S, Se. A ball-and-stick model of this structure is 
shown in Fig. [T] Besides its natural occurrence in mineral 
form, stibnite can be synthesised using a variety of low- 
cost fabrication techniques J£r— Using these techniques 
it is possible to obtain a good degree of crystallinity ) 19 i 20 
to control dimensionality ^ 16 ' 21 and to tune the optical 
properties. 16 ' 22 ~ — Semiconductors of the stibnite fam- 
ily have also been synthesized in various nanostructured 
forms. For example Refs. I5I H4IIT5I and Refs. I5II24I reported 
nanowires and nanotubes, respectively, of stibnite, anti- 
monselite and bismuthinite. Nanowires of stibnite were 
found to exhibit enhanced ferroelectric and piezoelcc- 




FIG. 1: Ball-and-stick model of A2B3 semiconducting metal 
chalcogenides of the stibnite family, with A standing for Sb or 
Bi (brown), and B for S or Se (yellow). The two inequivalent 
(A4Bs) n ribbons in the unit cell are highlighted in red, and 
the perspective view is along the direction of the ribbons. 
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trie properties as compared to their bulk counterpart j^ 5 . 
Nanowires and nanotubes of antimonselite were found 
to exhibit conductivities much higher than their bulk 
counterpart^ and are being considered for thermoelec- 
tric applications. In the case of bismuthinite, Ref. [l4| 
reported nanowires with diameters as small as 1.6 nm, 
corresponding to a transverse size of only two ribbons. 
The rhombohcdral phase of Bi 2 Sc 3 has been investigated 
extensively since this compound is a prototypical topo- 
logical insulator^ However to the best of our knowledge 
little is known about orthorhombic Bi 2 Se3, i.e. guana- 
juatite, which is stable only at high temperature and 
pressure. 27 ' 28 

The band gaps of stibnitc, antimonselite and bismuthi- 
nite have been measured extensively via optical absorp- 
tion experiments. The band gap of stibnite ran ges be- 
tween 1.42-1.78 eV, 2 ^ For antimonselite Ref.^H re- 
ported a direct gap of 1.55 eV, while Ref. HU gave an 
indirect gap between 1-1.2 eV. The measured band gap 
of bismuthinite is 1.38-1.58 eVM^- The spread in the 
measured gaps can be attributed to the different prepa- 
ration conditions used, yielding different degrees of poly- 
crystallinity and even amorphous samples in some cases, 
and also different stoichiometrics. In addition all these 
compounds exhibit closely lying direct and indirect tran- 
sitions (cf. Fig. [3] below) , thereby complicating the as- 
signment of the nature of the optical gap. 
All four minerals of the stibnite family have been investi- 
gated in detail using density-functional theory (DFT) cal- 
culations. The electronic properties of these compounds 
have been studied in Refs. ill-Ill and the elastic and 
optical properties have been calculated in Ref. [3^. 
A comparison of the theoretical studies published so far 
shows some inconsistencies in the calculated band gaps, 
for example the values reported for stibnite are in the 
range 1.18-1.55 eV&2&22^ As expected all the calcu- 
lated DFT gaps underestimate the measured band gaps. 
To the best of our knowledge only one work^ 3 . reported a 
calculation of the quasiparticle band gap of stibnite and 
antimonselite within the GW approximation^ 
Within this context there exists a need for detailed and 
reproducible calculations of the electronic structure of 
stibnite and related compounds based on state-of-the- 
art quasiparticle techniques. In line with this need the 
goal of the present work is to report a systematic and re- 
producible study of the quasiparticle band structures of 
all four A2B3 semiconducting metal chalcogenides of the 
stibnite family. An emphasis is placed on convergence 
tests and on the sensitivity of the quasiparticle correc- 
tions to the structural parameters and the inclusion of 
semicore d states in the calculations. 
Our calculated G0W0 band gaps are 1.6 eV (Sb 2 S 3 ), 1.3 
eV (Sb 2 Se 3 ), 1.7 eV (Bi 2 S 3 ) and 1.2 eV (Bi 2 Se 3 ). These 
values differ from the corresponding DFT gaps by up to 
0.5 eV, and are close to the measured optical gaps. By 
inspection of the band structures we infer that all four 
compounds have direct band gaps, although in most cases 
an indirect transition just below the direct gap (within 



0.1 eV) is also possible. The inclusion of semicore elec- 
trons in the calculations is found to modify the band 
gaps by 0.1-0.2 eV. In addition we find that the gaps are 
rather sensitive to the lattice parameters, as they change 
by up to 0.3 eV when the lattice parameters are taken 
from experiment or fully optimized within DFT. 

The manuscript is organized as follows. In Sec. [TT] we 
describe the computational methodology and the con- 
vergence tests for the GW calculations. In Sec. IIIII we 
present our main results, including quasiparticle band 
structures and band gaps. In Sec. IIVI we discuss our 
findings in relation to the photovoltaics applications of 
the materials considered in this work. In Sec. [V] we sum- 
marize our results and present our conclusions. 



II. METHODOLOGY 

A. DFT calculations 

All DFT calculations are performed using the Quantum 
Espresso packaged The calculations are based on the 
local density approximation (LDA) to DFTj 46 ' 47 
Only valence electrons are explicitly described, and the 
core- valence interaction is taken into account by means of 
Troullier-Martins pseudopotentials 4 ^ generated using the 
f hi98 code42 In the cases of S (Se) the 3s 2 3p 4 (4s 2 V) 
electrons are included in the valence as usual. For Sb 
and Bi we generate two sets of pseudopotcntials, one set 
with five electrons in the valence, i.e. 5s 2 5p 3 and 6s 2 6p 3 
respectively, and one set with additional semicore 4<i 10 
and 5d 10 electrons, respectively. 

The electronic wavefunctions are expanded in planewaves 
basis sets with kinetic energy cutoffs of 70 Ry (Sb 2 Se 3 , 
Bi 2 Se 3 ) and 90 Ry (Sb 2 S 3 , Bi 2 S 3 ) for the calculations 
without semicore states, and 100 Ry (Bi 2 S 3 , Bi 2 Se 3 ) and 
130 Ry (Sb 2 S 3 , Sb 2 Se 3 ) when semicore states are in- 
cluded. In each case considered the selected cutoff yields 
a total energy converged to within 2 meV/atom. All 
self-consistent calculations are carried out using a 8x8x8 
Brillouin zone mesh centered at T, corresponding to 170 
irreducible points. 



B. Crystal structure 

Stibnitc (Sb 2 S 3 ), antimonselite (Sb 2 Sc 3 ), bismuthinite 
(Bi 2 S 3 ), and guanajuatite (Bi 2 Se 3 ) all crystallize in the 
same orthorhombic lattice and belong to the Pnma 62 
space group^ Each unit cell contains 20 atoms, whose 
coordinates can be generated by applying the symme- 
try operations of the crystallographic group to a set of 
5 atomic coordinates. Figure Q] shows a ball-and-stick 
representation of these A 2 B 3 structures. The structural 
parameters were measured by Refs l27ll50l - [52l and are re- 
ported in Ref. HH. 
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C. Quasiparticle calculations 

We calculate the quasiparticle energies within many- 
body perturbation theory using the GW method ) 44 ' 53 " — 
as implemented in the SaX code£l The GW self-energy 
is evaluated in the GoWo approximation as E = iGqWo- 
Here Go denotes the electron Green's function defined by 
the Kohn-Sham eigenstates ip n \s.(r) and eigenvalues e„k 
corresponding to the band index n and the wavevector 
k, and Wo represents the screened Coulomb interaction 
calculated in the random phase approximation . 53 ' 58 The 
quasiparticle energies E n k are obtained as&& 

E n k = e„k + Z„ k (V>„k|E(e nk ) - T4<#nk), (1) 

where E n k is the quasiparticle energy, Z n v. the associated 
quasiparticle renormalization, and V xc is the exchange 
and correlation potential. 

The self-energy can be written as the sum of a bare ex- 
change contribution E x and a correlation contribution 
E c : E = E x + E c . The exchange part does not depend 
explicitly on the excitation energy and readsi 5 -^ 

E x (r,r') = -£ Ck(r)V>„ k (i>(r,r'), (2) 

n£occ,k 

where the sum is over occupied states and v represents 
the bare Coulomb interaction. This contribution to the 
quasiparticle correction is sensitive to the overlap be- 
tween Kohn-Sham wavefunctions regardless of their en- 
ergy. As a result the use of semicore states can have sig- 
nificant effect on the calculations, as shown in Refs. |60| - 
|H. This aspect will be discussed in detail in Sec. IHI CI 
The energy-dependence of the correlation contribution 
E c arising from the dynamically-screened Coulomb in- 
teraction is described via the Godby-Needs plasmon-pole 
model.— We use a plasmon-pole energy of 1 Ry for all 
materials, similar to the energy of the peaks in the cal- 
culated electron energy loss spectra. 
Since the computational efforts for achieving convergence 
in E c and E^ are very different owing to the necessity of 
evaluating unoccupied states for E c , we perform sepa- 
rate convergence tests for these two components. For 
the exchange contribution we use kinetic energy cutoffs 
of 75 Ry and 100 Ry for calculations without and with 
semicore electrons, respectively. For the correlation con- 
tribution we perform convergence tests by calculating the 
band gap at various kinetic energy cutoffs up to 7 Ry 
for the polarizability. Figure HJa) shows that the band 
gap is converged within 0.05 eV already for a cutoff of 
5 Ry. The dependence of the band gap on the polarizabil- 
ity cutoff shows the same trend for calculations with or 
without semicore states. This is consistent with the ex- 
pectation that the effect of semicore states in E c should 
be smalls Based on the data of Fig. [U[a), in the fol- 
lowing we use a polarizability cutoff of 7 Ry for calcu- 
lations without semicore electrons, and of 6 Ry for the 
more demanding calculations including semicore states. 
In Fig. HJb) we show the convergence of the band gap of 



antimonselite with respect to the energy of the highest 
unoccupied state included in the polarizability. Based 
on the trend in this figure we set the number of unoc- 
cupied states to 224 and 264 for calculations with and 
without semicore, corresponding to a maximum energy 
denominator of 35 eV. Both Go and Wq are calculated on 
uniform and T-centered 2x6x2 Brillouin-zone meshes. 
In order to estimate the accuracy of our quasiparticle cor- 
rections with respect to the above convergence parame- 
ters we follow the approach of Ref. H3- In this approach 
the dependence of the band gap on a given convergence 
parameter is fitted by a simple function in order to ex- 
tract a "best-guess" asymptotic limit. This asymptotic 
limit is then taken to represent the converged gap. In this 
work we tentatively approximate gap vs. cutoff curves us- 
ing the following function: 

#g Q ap = a +a 1 (x-a 2 )- 1 / a \ (3) 

where E^ is the quasiparticle band gap, x is the con- 
vergence parameter (i.e. the polarizability cutoff or the 
largest energy denominator) and clq, . . . ,0,3 arc fitting pa- 
rameters. While Eq. (|3|) is largely arbitrary, this choice 
reflects the expectation that the gap will converge faster 
than 1 jx owing to the damping introduced by the matrix 
elements in the Adler- Wiser polarizability 65 ' 66 Figure [H 
shows that the fitting curves obtained for stibnite de- 
scribe rather accurately the calculated data points, there- 
fore it is reasonable to assume that the parameter a>o ob- 
tained from the fit should provide a good estimate of the 
converged gap. By repeating this procedure for all four 
compounds Sb2S3, Sb2Se3, Bi2S3, and Bi2Se3 we find 
that the convergence parameters described above yield 
band gaps which differ by less than 0.05 eV from the 
corresponding asymptotic values. 



III. RESULTS 

A. Structural parameters 

Table U shows the comparison between our calculated lat- 
tice parameters and experiment. As expected the use of 
the local density approximation leads to a general under- 
estimation of the experimental parameters. Interestingly, 
while such underestimation does not exceed 1.1% along 
the direction of the (A4Be) n ribbons (b parameter in Ta- 
ble |U cf. Fig. [TJ, the deviation can reach up to 4.2% in 
the direction perpendicular to the ribbons (a and c pa- 
rameters in Table H]). We tentatively assign this behavior 
to the fact that inter-ribbon forces are likely to include 
non-negligible van der Waals components, and hence are 
not described correctly within the LDA. Inspection of the 
calculated cohesive energies seems to support this possi- 
bility. In fact, while we calculate an intra-ribbon cohe- 
sive energy of 4.1-4.8 eV/atom for the four compounds 
considered, the inter-ribbon cohesive energy is only of 
0.2-0.3 eV/atom. 
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FIG. 2: (color online) Convergence tests for the quasiparticle band gap of antimonselite. (a) Calculated Go Wo band gap as a 
function of the polarizability cutoff, for calculations without (red circles) or with (blue disks) semicore states. The solid lines 
correspond to the fits obtained from Eq. @- We find ao = 1.33 eV and 1.16 eV for calculations without and with semicore 
electrons, respectively, (b) Same as in (a), with the band gap reported as a function of the largest energy denominator used 
for calculating the polarizability. In this case we find oo=1.32/1.14 eV for calculations without/with semicore electrons. All 
calculations were performed using optimized lattice parameters. 



In order to take into account these differences between 
calculated and experimental lattice parameters, in the 
following we report quasiparticle calculations obtained 
for both sets of parameters. 



B. DFT/LDA band structures 

Figure [3] shows the DFT/LDA band structures calculated 
using experimental lattice parameters and without semi- 
core electrons. Calculations including the semicore states 
yield very similar band structures. For clarity we only 
show the dispersions along the Z-T-X path and along 
the Y-T segment running along the axis of the (A4Be)„ 
ribbons. The top of the valence band is found to be 
predominantly of S-3p or Se-4p character, while the bot- 
tom of the conduction band comprises of Sb-5p or Bi-6p 
states, consistently with previous calculations ! 36 ' 42 

The band structures shown in Fig. [3] exhibit several ex- 
trema in proximity of the fundamental gap, making the 
direct and indirect transitions almost degenerate. Ta- 
ble |ll] shows that the energy separation between direct 
and indirect DFT/LDA band gaps falls within 0.15 eV 
in all cases. The data in the table suggest that in these 
compounds the direct transition will most likely dom- 
inate over the indirect one, apart from a very narrow 
onset of 0.1-0.2 eV. This observation is consistent with 
experimental evidence showing a weak absorption onset 
just below the threshold for direct absorption ,11*211 There- 
fore for practical purposes, and in particular for photo- 
voltaics applications, stibnite, antimonselite, bismuthi- 
nite and guanajuatite can be considered as "effectively 
direct gap" semiconductors. 



C. Quasiparticle corrections 

Figure [3] shows that GW quasiparticle corrections lead 
to a moderate increase of the band gaps in all cases, 
while generally preserving the shape of the band extrema. 
From this figure we deduce that a simple scissor operator 
should be able to capture the most important effects of 
the GW corrections. 

A detailed analysis of the quasiparticle corrections at the 
high-symmetry points T, X, and Z is given in Fig. 0] 
and Table IIII1 In Fig. [4] we report the quasiparticle cor- 
rections as a function of the corresponding Kohn-Sham 
eigenvalues around the band extrema. In the cases of 
stibnite and antimonselite we observe that in the cal- 
culations with semicore electrons the valence bands are 
slightly upshifted (by about 0.1 eV) as compared to cal- 
culations without semicore, while the corrections to the 
conduction bands are essentially the same. In the cases 
of bismuthinitc and guanajuatite the effect of semicore is 
to shift the valence bands up and the conduction bands 
down by a similar amount (~0.1 eV). As a result of these 
small changes, the quasiparticle corrections to the band 
gaps calculated with or without semicore electrons can 
differ by up to 0.2 eV (cf. Table HIIJ). 
Semicore electrons appear to slightly reduce the quasi- 
particle corrections as compared to calculations without 
the semicore. This finding is consistent with previous 
calculations and can be rationalized as follows^2r— The 
semicore d states introduce additional contributions 
and c to the GW self energy. Of these contributions, 
while the correlation part is small owing to the large 
energy separation between semicore states and conduc- 
tion states, the exchange part can be large since it 
does not contain energy denominators but is sensitive to 
the overlap between the band edge states and the semi- 
core states. This interpretation is confirmed by Fig. (SJ 
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Experiment 



Calc. w/o semicore Calc. with semicore 



Sb 2 S 3 11.311° 3.836° 11.229° 11.036 3.795 10.753 11.087 3.838 10.834 

-2.4% -1.1% -4.2% -2.0% 0.1% -3.5% 

Sb 2 Se 3 11.794° 3.986° 11.648'' 11.609 3.952 11.213 11.646 3.989 11.287 

-1.6% -0.9% -3.7% -1.3% 0.1% -3.1% 

Bi 2 S 3 11.305° 3.981° 11.147° 11.227 3.999 11.001 11.030 3.949 10.853 

-0.7% 0.5% -1.3% -2.4% -0.8% -2.6% 

Bi 2 Se 3 11.830° 4.090 d 11.620 d 11.767 4.141 11.491 11.609 4.099 11.374 

-0.5% 1.3% -1.1% -1.9% 0.2% -2.1% 

° Ref. HO- 
b Ref. Hl|. 
° Ref. |H. 
d Ref. H3. 

TABLE I: Comparison between the calculated DFT/LDA lattice parameters of stibnite, antimonselite, bismuthinite, and 
guanajuatite and experiment (all values are given in A). The percentile deviation from experiment is indicated in each case. 
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FIG. 3: (color online) Band structures of (a) stibnite, (b) antimonselite, (c) bismuthinite, and (d) guanajuatite calculated using 
DFT/LDA, experimental lattice parameters, and without semicore electrons (black solid lines), as well as corresponding density 
of states (DOS, black dashed lines). The contributions to the DOS from the p states of S and Se (Sb and Bi) are indicated 
by the green (blue) shaded areas in each case. The GW quasiparticle energies of the band extrema at high symmetry points 
are also shown, with blue squares and red circles indicating calculations with or without semicore electrons, respectively. The 
connecting lines are guides to the eye. The coordinates of the high symmetry points in reciprocal lattice units are as follows: 
Z : (0,0,0.5), X : (0.5,0,0), Y : (0,0.5,0). 
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FIG. 4: Quasiparticle corrections as a function of the corresponding DFT/LDA eigenvalues for (a) stibnite, (b) antimonselite, 
(c) bismuthinite, and (d) guanajuatite. Only eigenvalues at the high-symmetry points Y, X and Z are considered. Blue disks 
and red circles indicate calculations with and without semicore electrons, respectively. 



Mimimum gap Direct gap 

w semicore w/o semicore w semicore w/o semicore 



Sb 2 S 3 


1.19 


1.21 


1.26 


1.27 


Sb 3 Se 3 


0.84 


0.86 


0.84 


0.86 


Bi 2 S 3 


1.25 


1.24 


1.28 


1.27 


Bi 2 Se 3 


0.85 


0.86 


0.99 


0.99 



TABLE II: Comparison between the minimum band gaps and 
the direct band gaps of stibnite, antimonselite, bismuthinite, 
and guanajuatite, as obtained from DFT/LDA. In these cal- 
culations we use the experimental lattice parameters. All val- 
ues are in units of eV. 



where we can see that the inclusion of semicore electrons 
does indeed affect the exchange part of the GW correc- 
tions, while at the same time the correlation component 
remains almost unchanged. 

Generally GW quasiparticle corrections result in nega- 



tive shifts of the valence band DFT/LDA eigenvalues^ 
At variance with this expectation in the systems con- 
sidered here we find that the corrections to the valence 
are mostly positive (cf. Fig. 2]). Cross-checking our re- 
sults in a different computational framework yields the 
same trends This unusual trend can be linked with the 
remarks made for CdS in Ref. |60| . In that work the au- 
thors compared the GW quasiparticle energies of CdS 
obtained for calculations using Cd pseudopotentials in- 
cluding a progressively higher number of semicore states. 
Their tests showed that the separation of the n = 4 shell 
into core s and p states and semicore d states led to 
significant inaccuracies in the exchange part of the self- 
energy corrections, resulting in unusual upshifts of the 
Cd-Ad states. This problem was resolved by including 
the entire n — 4 shell in the semicore. In the present 
work we only include the d states of the n = 4 shell of 
Sb and the n = 5 shell of Bi, therefore a small system- 
atic inaccuracy in the absolute value of the quasiparticle 
shifts cannot be excluded. We note however that in our 
case the states near the band edges are of p-character, 
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FIG. 5: Exchange (black disks and circles) and correlation (red filled and empty squares) contributions to the quasiparticle 
corrections vs. DFT/LDA eigenvalues for (a) stibnite, (b) antimonselite, (c) bismuthinite, and (d) guanajuatite. Only eigenval- 
ues at the high-symmetry points F, X and Z are considered. Filled and empty symbols indicate calculations with and without 
semicore states, respectively. 



therefore they should be less sensitive to the mechanism 
identified in Ref. HO for the valence d states of Cd. We 
also point out that our calculations do not include the 
semicore states of Se, and this may also contribute to the 
upshift. Owing to the large number of atoms (20) in the 
unit cells of these compounds, testing these speculations 
would be highly impractical. 

Table Hill reports the DFT/LDA eigenvalues and the cor- 
responding quasiparticle corrections for the valence band 
top and conduction band bottom at the high symmetry 
points r, X and Z. From this table we see that the LDA 
band gaps at these points are sensitive to the choice of 
the lattice parameters, and this sensitivity is reflected 
in the corresponding quasiparticle energies. Calculations 
performed using optimized lattice parameters or experi- 
mental parameters can differ by up to 0.2 eV. This obser- 
vation may explain the lack of consensus between previ- 
ous computational investigations of the band structures 
of these compounds^ ' 36 ' 39 ' 43 

Taken together the sensitivity of the quasiparticle en- 
ergies to the presence of semicore electrons and to the 



choice of lattice parameters leads to differences in the 
calculated values which can be as large as 0.4 eV. This 
finding suggests that it is important to use extreme care 
when comparing the quasiparticle band structures of stib- 
nite and related compounds with experimental data. 
In the remainder of this manuscript we will focus on cal- 
culations using experimental lattice parameters and in- 
cluding semicore electrons, which we consider our best 
estimates for the quasiparticle energies in these com- 
pounds. 



D. Band gaps 

Table IIVI reports the band gaps calculated using the 
experimental lattice parameters and including semicore 
electrons. The band gaps are obtained by considering 
the band extrema at T, X, and Z , and we give both the 
fundamental gap and the direct gap. While in stibnite, 
antimonselite, and bismuthinite the calculated minimum 
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GW gap is indirect, the difference between direct and 
indirect gaps is within 0.1 eV. In guanajuatite the fun- 
damental GW gap is direct. These results suggest that 
all four compounds can be considered direct-gap semicon- 
ductors for practical applications, especially in the area of 
optoelectronics. The calculated direct gaps are 1.60 eV 
(stibnite), 1.28 eV (antimonselite) , 1.72 eV (bismuthi- 
nite), and 1.18 eV (guanajuatite). As shown in Tabic ITV 
these values are in line with previous GW calculations 
where available ; 37 ' 43 and also rather close to measured 
optical gaps. 

The comparison with experimental data is not straight- 
forward since the experimental literature appears to only 
report optical gaps (cf. literature review in Tabic liV|) . 
However our calculations refer to quasiparticle gaps and 
do not include excitonic effects. Including excitonic 
effects using the Bcthc-Salpctcr approach^ would be 
rather challenging owing to the large size of these sys- 
tems. To the best of our knowledge no excitonic effects 
were measured or mentioned for any of the four com- 
pounds studied. One exception is possibly the absorp- 
tion spectrum reported in Ref. [Til which exhibits some 
sharp features resembling excitonic peaks, however the 
authors assigned those peaks to defects or internal reflec- 
tions. The agreement between our calculated quasipar- 
ticle gaps and the measured optical gaps can be seen a 
posteriori as an indication that excitonic shifts are small 
in this class of semiconductors. 

Figure [5] provides a schematic view of our calculated 
quasiparticle band gaps, and shows a comparison with 
DFT/LDA calculations and experiment. 

IV. DISCUSSION 

Taking the calculated quasiparticle band gaps of 1.2- 
1.7 eV as representative of the optical gaps, the four 
semiconductors considered here lie precisely in the range 
of the optimal Shocklcy-Queisscr performance^ The 
Shockley-Queisscr analysis addresses the ultimate effi- 
ciency of a solar cell based on a single material as light 
absorber and electron conductor, e.g. silicon solar cells. 
In this analysis the optimum efficiency results from a 
trade-off between maximizing the band gap in order to 
increase the photovoltage, and minimizing the band gap 
in order to increase the photocurrenti^ 
In the case of nanostructured solar cells based on 
the donor/acceptor concept such as for instance 
semiconductor-sensitized solar cells; 1 -' 7 - the Shocklcy- 
Queisser analysis needs to be modified in order to 
take into account the energy-level alignment at the 
donor/acceptor interface. In fact, while the photocur- 
rent is still determined by the optical gap of the ab- 
sorber (typically the donor), at variance with conven- 
tional bulk solar cells the photovoltage is dictated by the 
difference between the lowest unoccupied states of the 
acceptor and the highest occupied states of the donor. 
This effect can be taken into account by introducing the 



Optimized parameters Expt. parameters 

LDA GW LDA GW 

w/o S S w/o S S w/o S S w/o S S 



TABLE III: Quasiparticle energies of stibnite, antimonselite, 
bismuthinite and guanajuatite at the high-symmetry points 
F, X, Z vs. the corresponding DFT/LDA eigenvalues. We 
report both sets of results obtained using optimized or ex- 
perimental lattice parameters. The columns labelled "S" and 
"w/o S" indicate calculations with and without semicore elec- 
trons, respectively. For each high-symmetry point we consider 
the energies at the valence band top (e.g. r„) and the con- 
duction band bottom (e.g. T c ). All values are in units of eV. 



concept of "loss-in-potential" ^ which is the reduction 
of the photovoltage resulting from the energy mismatch 
and additional losses. Loss-in-potentials estimated for 
actual devices can be as large as ~ 1 eV, and the most 
optimistic scenario would correspond to losses as small 
as 0.3 eVi 7 ^ Figure [7] shows the theoretical efficiency of 
semiconductor-sensitized solar cells based on stibnite, an- 
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0.10 0.20 0.00 

1.58 1.52 1.33 
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1.65 1.65 1.22 
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0.07 0.18 -0.11 

1.40 1.32 1.09 

0.19 0.29 0.00 

1.53 1.46 1.10 
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1.37 1.35 0.94 

Bismuthinite 

-0.04 0.08 -0.14 

1.57 1.48 1.28 

0.02 0.09 0.00 

2.04 2.09 1.61 

-0.13 -0.04 -0.08 

1.89 1.93 1.41 

Guanajuatite 

0.18 0.30 -0.07 

1.39 1.24 1.17 

0.19 0.28 0.00 

1.61 1.58 1.04 

-0.06 0.05 -0.14 

1.63 1.65 1.16 
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1.42 (1.30) 




1.72 (1.64) 


1.4°, 1.38 p , 1.58 g ' j 


Bi 2 Se 3 


0.9 a , l.l r 


1.06 




1.18 





a Rcf.Hi, 6 Rcf.lH c Ref.H d Ref.lH, e Ref.ljl / Ref.lTJ 

9 Ref.lH, h Rcf.113, i Ref. El j Ref.[Il fe Ref. El ; Ref . [3l|, m Ref.ll 

n Ref. lU ° Ref. lH, " Ref. lH, i Ref. El r Ref. |4(] 



TABLE IV: Comparison between calculated and measured band gaps of stibnite, antimonselite, bismuthinite, and guanajuatite. 
We report the direct band gaps calculated within DFT (Kohn-Sham gaps) and GW (quasiparticle gaps), and the measured 
optical gaps. The values in parenthesis indicate the calculated indirect band gaps in each case. All band gaps are in eV. The 
direct gaps are reported for Y (Sb 2 S 3 and Bi 2 S 3 ), X (Bi 2 Se 3 ), and Z (Sb 2 Se 3 ). Our calculations include semicore electrons and 
were performed using the experimental lattice parameters. 



timonselitc, bismuthinite, and guanajuatite, calculated 
using the prescription of Ref. [70| for a loss-in-potential 
of 0.3 eV. While these estimates are very crude and the 
projections are possibly too optimistic, it is interesting 
to note that all of these four materials cluster very near 
the optimum power conversion efficiency of 20-25%. 
From Fig. [7] we infer that the four compounds studied 
here are all promising candidate for nanostructured pho- 
tovoltaic applications, with antimonselite and guanaju- 
atite slightly superior to stibnite. In particular it can- 
not be excluded that guanajuatite, even if unstable at 
room temperature in bulk form, could be stabilized as a 
nanostructure. Given its projected maximum efficiency 
in Fig. [3 it might be worth to attempt the synthesis of 
guanajuatite nanoparticles. In the case of bismuthinite 
Refs. [fUH showed that this material does not work as a 
semiconductor sensitizer for Ti02, owing to the incor- 
rect energy-level alignment at the interface. However it 
cannot be excluded that bismuthinite could still reach 
the ideal efficiency when combined with an alternative 
acceptor, e.g. SnC>2 or ZnO. 



V. CONCLUSIONS 

In this work we report a systematic study of the quasi- 
particle band structures of the four isostructural metal 
chalcogenides stibnite (Sb2S3), antimonselite (Sb2Se3), 
bismuthinite (Bi2S3), and guanajuatite (Bi2Se3), within 
the GW approximation. 

In order to ensure reproducibility of our results we have 
placed an emphasis on convergence tests and explored the 
effects of various calculation parameters, such as for in- 
stance the role of semicore d electrons and lattice param- 
eters. The inclusion of semicore electrons in the calcula- 
tions is found to modify the band gaps by up to 0.2 eV, 
and the choice of experimental vs. optimized lattice pa- 
rameters can lead to differences of up to 0.2 eV in the cal- 
culated gaps. These findings indicate that some caution 
should be used in discussing the theoretical band gaps of 
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SbfeS, Sh.Se, Bi& Bi,Se, 



FIG. 6: Schematic summary of the band gaps of stibnite, 
antimonselite, bismuthinite, and guanajuatite calculated in 
this work: Kohn-Sham gaps (empty rectangles) and GW gaps 
(filled rectangles). The band gaps were obtained by including 
semicore electrons and using the experimental lattice param- 
eters. The pink rectangles indicate the range of experimental 
optical gaps reported in Table [TVl 
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Band gap (eV) 

FIG. 7: Ideal efficiency of nanostructured solar cells based 
on semiconductors of the stibnite family. The theoretical ef- 
ficiency as a function of the band gap energy (black solid 
curve) is calculated using the prescription of Ref. [7^ with a 
loss-in-potential of 0.3 eV and a fill factor of 73%. 
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these materials and in comparing with experiment. 
Our calculations indicate that all four compounds have 
direct band gaps, barring indirect transitions marginally 
below the direct gap. The calculated band gaps are 
1.6 eV (stibnite), 1.3 eV (antimonselite), 1.7 eV (bis- 
muthinite) and 1.2 eV (guanajuatite). These values fall 
within the range of measured optical gaps, although it 
must be observed that there is a considerable scatter in 
the experimental data, possibly due to different prepara- 
tion conditions. 

Using a modified Shockley-Queisser analysis^ we es- 
timate the ultimate performance of solar cells based 
on these compounds as light sensitizers. This analy- 
sis indicates that all four materials have potential for 
high-efficiency nanostructured solar cells. The highest 
theoretical efficiency is obtained for guanajuatite, the 
high-temperature polymorph of the topological insula- 
tor Bi2Se3. This finding suggests that it could be worth 
attempting to synthesize nanoparticles of this compound 
(which is unstable in bulk form at room temperature). 
Future calculations should address the optical absorption 



spectra of these compounds within the Bcthc-Salpctcr 
approach, in order to establish whether excitonic effects 
are small as our data appear to suggest. It will be also 
interesting to extend the present study to the case of in- 
dividual nanoribbons of these metal chalcogenides, since 
liquid-phase exfoliation techniques for van der Waals 
bonded materials are becoming increasingly popular^ 

We hope that the present study will contribute to the on- 
going research on new materials for energy applications, 
and stimulate further efforts to understand and exploit 
these fascinating and relatively unexplored compounds. 
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